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Abstract. A monolayer graphene exists in an environment where a uniform magnetic 
(~| , field interacts a spatially modulated magnetic field. The spatially modulated magnetic 

^ ' field could affect Landau levels due to a uniform magnetic field. The modulation 

D . effects on Landau levels are investigated through the Peierl's tight-binding model. 

The magneto-electronic properties are dominated by the period, the strength, and the 
direction of a spatially modulated magnetic field. Such a field could induce the growth 
in dimensionality, the change of energy dispersions, the destroy of state degeneracy, 
and the creation of band-edge states. There are a robust Landau level at Fermi level 
. and ID parabolic subbands located around the original Landau levels, which make 

C I density of states exhibit a delta-function-like structure and many pairs of asymmetric 

peak structure, respectively. The density of states and the energies of band-edge states 
strongly depend on the strength, but not on the period and the direction. 
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Recently, the discovery of few- layer graphene [H [2] has inspired a lot of theoretical 
and experimental studies in condensed matter physics and material science. Graphene, 
as a nanomaterial, not only has its academic value in condensed matter physics, but 
also promises to be the candidate for the next-generation transistor [1]. Therefore, 
any advances in understanding its elementary physical properties, such as electronic 
properties [SliaElElITllHlElIO], transport properties [IIl[l2l[l3l[Tl[l5l[l6l[I71[l8l[l9] 
, optical properties [20l [211 |22l [231 [2l] and many body effects |l|25l [26] would be very 
valuable for finding possible applications. The fact that few-layer graphene is accessible 
from experiments means that they can be controlled in single-atom thick accuracy. From 
a technical point of view, this motivates the idea of engineering the electronic structures 
[H [3] . The main reason is that the electronic properties of few-layer graphene strongly 
depend on the number of layers and the stacking order [6l H]. One can further tune 
its electronic properties by applying a gate voltage. Such an external field has been 
shown to induce the transition of fundamental carriers between electrons and holes in a 
monolayer graphene [1] , and lead to opening an energy gap in a bilayer graphene [31 [7] . 

The fact that Graphene continuously surprises us should owe to its exotic electronic 
structure, where low-energy quasiparicles behave as massless Dirac particles [H]. In a 
uniformly perpendicular magnetic field, the electronic states flock together and form 
unusual OD Landau levels. The level energies, unlike those of 2D electron g as , are 
proportional to the square root of the quantum number and the field strength, rather 
than being equally spaced or proportional to the field strength. Recently, these features 
have been verified by measurements based on the infrared spectroscopy [20] . The Landau 
level quantization reflects the Dirac nature of its quasiparticle, and is deduced to be 
responsible for the observation of unconventional integer quantum Hall effect [HI [13]. 
On the other hand, the spatially modulated electric and magnetic fields were predicted 
to cause drastic changes in state degeneracy, energy dispersion, band-edge state, and 
band width [271 EH]- In this work, we would like to study the modulation effects on 
Landau levels due to a spatially modulated magnetic field. The influences of the field 
strength, period, and direction on Landau levels will be investigated in detail through 
employing Peierl's tight-binding model. 

Graphene is a honeycomb lattice with hexagonal symmetry. The vr-band structure, 
formed by 2pz orbital from each carbon atom, can be calculated through the tight- 
binding model within the nearest-neighbor atomic interactions. Since there are two 
carbon in the primitive cell, the Hamiltonian can be represented by a 2 x 2 Hermitian 
matrix in the space expanded by Bloch functions of two crystalline sublattices. To apply 
tight-binding formalism to Bloch electrons in a magnetic field, we need to consider 
the extra phase factor (Peierl's phase), which depends on the vector potential A, in 
describing the wave function. Within the tight-binding scheme [29]. the wave function 
in a magnetic field is expressed as 



R„ is the lattice vector with sublattice index n, ipp^ is the atomic orbital, and 




(1) 
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= A(^) ■ d$, is the phase associated with the magnetic field. It is obvious 
that, by disregarding Gr„, the wave function can recover the original Bloch function. 
In fact, the phase {e/fijCn would lead to the breakdown of the Bloch condition for the 
wave function. With this modified Bloch function the Hamiltonian matrix element 
has already been shown to be that in zero field multiplied by a phase factor in such a 
condition that the magnetic field changes slowly as a function of the lattice constant. 
It is given by 

^kT' = (*kl^B|0 

= ^ 5] exp{-z(k • R„ - k' • R„0 - ^^i'^n, Rn')}^RnR„M (2) 

where i^R^R^, = (y9p^(r— R„)|p^/2m+y|(/9p^(r— R„/)) is the Hamiltonian matrix element 
for a single particle in a periodic lattice potential V, and is nonvanishing for = 
R„+a„ (a„ is the position vector that connects an atom at R^ and its nearest neighbors). 
The Hamiltonian of Bloch electrons in a magnetic field is different from that in zero 
field, in the existence of the phase difference term Q(R„,R„/) = {e/h){Gn^ — Cr^,). 
J1(R„, R„/) could give rise to the couplings between states with different k's, and hence 
make solving band structure more complicated. Generally, the Hamiltonian in k space 
would be irreducible. In other words, the magnetic field completely destroys the crystal 
symmetry and k is not a good quantum number anymore. However, if r2(Rn, Rn') is 
periodic function of R„, the periodicity of crystal will be preserved, and the dimension 
of the unit cell is in turn determined by the period of fl. 

The graphene exists in an environment where a uniform perpendicular magnetic 
field Bz interacts a spatially modulated magnetic field B' sin Kxz along the armchair 
direction with the period Ib — 27: /K [Figure 1]. The magnetic fiux, the product 
of magnetic field and the hexagonal area, is $ = 3y/3b^/2 in unit of fiux quantum 
($0 = hc/e = 4.1356 x 10^^ [T/m^]). b = 1.42A is the C-C bond length. The vector 
potential can be chosen as A = {Bx — B' cos Kx/ K)y. Such a vector potential would 
lead to the new periodicity along the armchair direction. The unit cell is thus enlarged 
and its dimension is determined by $ = 1/Rb and Ib = 3bR'^. The parameters 
Rb and R'^ are both chosen as positive integers to accommodate the requirement of 
finite dimensionality of the Hamiltonian matrix, or finite dimension of a unit cell. The 
dimensionality is determined by the least common multiple {2Rm) of 2Rb and 2R'^. 
Therefore, the rectangular unit cell contains 4it!^ carbon atoms and the Hamiltonian 
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matrix is a 4i?m x 4:Rm Hermitian matrix, given by 
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Pn = iik(^) +hk{n) and q = igk- tik(^) = 7o exp[ (iA;^.6/2 + iky^/3b/2) + G^], t2k(^) = 
7o exp[ (iA;2,6/2 — iky\^b/2) — and tsk = 7oexp( — ifc^-fe) are three nearest-neighbor 
atom-atom interactions, where G„ = [in^^n — 1) + 1/6] + [—i[6{R'^y^'/n] cos[7r(n — 
5/6)/-Ry sin(7r/6-R^) ] is the phase caused by the magnetic fields. The magnetic flux 
$' = 3V3B'b'^/2, due to the modulated magnetic field, is used to characterize its 
strength. Note that we have adapted special arrangement of base functions to represent 
the Hamiltonian matrix as a band-like matrix. The base functions in the unit cell are 
chosen from outside to inside rather than from left to right [Figure 1] . The vr-electronic 
structure could depend on the direction of the modulated magnetic field, mainly owing 
to the anisotropy of graphene. If the magnetic field is modulated along the zigzag 
direction, the corresponding Hamiltonian matrix can be obtained in a similar way (not 
shown) . 

The unoccupied conduction bands (i^'^'s) are symmetric to the occupied valence 
bands [E'^^s) about the Fermi level Ep = 0. Only the former are discussed in this work. 
Because the range of is much smaller than that of ky for large Rm, it is sufficient 
just to consider the dispersion along ky in the following discussion. We first discuss 
the effects due to the modulated magnetic field along the armchair direction on the 
low energy bands resulted from the uniform magnetic field B = 40T. Without B', B 
could make electronic states fiock together and preform Landau levels E'^{n) (n is a 
nonnegative integer), as shown in Figure 2(a) by the dotted curves. Such levels are 
dispersionless and fourfold degenerate. These well-separated energy levels suggest that 
graphene under a uniform magnetic field could be regarded as a zero-dimensional system. 
In addition, the energies of Landau levels obey the simple relation E^{n) oc V uB. The 
perturbed modulated magnetic field {B' ^ 0) leads to the drastic changes in energy 
dispersions, degeneracies and band-edge states, as shown in Figure 2(a) by the solid 
curves at B' = AT with R = Rb'/Rb = 1- The dispersionless Landau levels with n > 1 
are changed into ID parabolic subbands, while the Landau level E^{n = 0) at Fermi 
level remains unaltered. Furthermore, such a field destroys the fourfold degeneracy and 
creates band-edge states. The magneto-electronic structure at high- and low-energy 
regimes are quite similar [Figure 2(b)]. At B' = 0, the high-energy Landau levels with 
E^{n) > 7o are doubly degenerate, and they are much closer to each other [30] than 
the low-energy Landau levels. Under the field [B' = 40 T, i? = 1), these Landau levels 
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are totally transformed into pairs of ID parabolic subbands. Each pair of parabolic 
subbands is nondegenerate and owns two band-edge states at the zone boundary. 

The strength, period, and direction of the modulated magnetic field strongly affect 
the low-energy bands, as shown in Figures 2(c)-2(f). The strength would change the 
band curvature and cause the shift in the energies of band-edge states. The energies of 
local maxima (minima) increase (decrease) in the increase of B', while the number (ric) 
of band-edge states keeps unchanged [Figure 2(c)]. The period presents diverse effects 
on the properties of the band structure, as shown in Figure 2(d) with shorter period 
and Figure 2(e) with longer period compared with Figure 2(a). As the period varies, 
the band curvature decreases when R < 1, whereas it is unaltered when R > 1. The 
number of subbands, however, shows different tendencies in the two situations. When 
R> 1, the number of subbands is proportional to R, while it is unchanged when R < 1. 
Despite the unchanged number of subbands, more band-edge states are created with 
He inversely proportional to R when R < 1. On the contrary, Uc is unchanged even 
though there are more subbands when R > 1. Furthermore, the period also affects 
state degeneracy in a particular manner. The energy bands are nondegenerate for most 
R's, but they are doubly degenerate as l/i? is an even integer. To see whether the 
anisotropy is of equal importance to the 7r-electronic structure, the modulated magnetic 
field with {B' = 4T,R = 1) is applied along the zigzag structure [Figure 2(f)]. Notice 
that the actual period in this case is VsERb' rather than SbRs' in armchair case. To 
make the comparison more precisely, one shall choose a proper Rb' of zigzag case to 
match the same period. Even with such consideration, it turns out that the influences 
of the modulated field along two definite directions on the Landau levels are similar to 
each other. 

Density of states (DOS), which directly reflects the main features of electronic 
structures, is defined as 



r(= 10~'^7o) is the phenomenological broadening parameter. The integration on k^, 
basically, could be neglected because of the very small range of k^. Without B', D{uo) is 
finite at Ep = and exhibits a lot of delta-function-like peaks at c<j 7^ 0, as shown by the 
dashed curve in Figure 3(a). Such symmetric prominent peaks come from the OD Landau 
levels at B=40 T. The distribution of peaks is nonuniform because of the unequally 
spaced Landau levels. When a spatially modulated magnetic field of {B' — AT, R — 1) 
is applied along armchair structure, every symmetric peak at a; 7^ is changed into a 
pair of asymmetric prominent peaks except for n = Landau level at Ep = 0, as shown 
by the solid curve in Figure 3(a). Such divergent structures come from band-edge states 
(maxima and minima) of ID parabolic subbands. Each pair of asymmetric prominent 
peaks is located around the energy of the original Landau level. The frequencies of the 
asymmetric prominent peaks are sensitive to the change of the strength [Figure 3(b)], 
whereas is insensitive to that of the period or the direction [Figure 3(c)]. The peak 
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height is closely related to the band curvatures about band-edge states and the number 
of band-edge states, both of which rely on the field strength and the period as mentioned 
earlier. The field strength could decrease the band curvature, so that the peak height 
gets weaker when B' increases. On the other hand, the period could decrease the band 
curvature and increase the number of band-edge states, in which the former makes the 
peak weaker, while the latter makes the peak stronger. The increasing compensates 
the decreasing, therefore, there is no net effect on the peak height. It is thus deduced 
that the low-energy spectrum of DOS is only affected by the strength of the modulated 
magnetic field even the period has complicated effects on the Landau levels. 

The frequencies of prominent peaks in DOS deserve a closer investigation. Figure 
4(a) shows the relation between the frequencies (cUc's) of the first six subpcaks and the 
strength at -B = 40 T and R= 1. The frequency of the first subpcak is fixed because its 
robustness against the modulated magnetic field as mentioned earlier. For other peaks, 
as the strength increases, the frequencies increase or decrease depending on whether the 
corresponding band-edge states are band maxima or minima. Similarly, the influence of 
the period is illustrated by its relation with (VcS of the first six subpeaks in Figure 4(b). 
The frequencies of those subpeaks are almost unchanged when the period varies. 

In summary, the modulation effects on Landau levels of a monolayer graphene 
are investigated through employing the Pcicrl's tight-binding model. The magneto- 
electronic properties are dominated by the strength, period, and direction of a spatially 
modulated magnetic field. Such a field could induce the growth in dimensionality, 
the change of energy dispersions, the destroy of state degeneracy, and the creation 
of band-edge states. The Landau levels are transformed into ID parabolic subbands, 
but a robust fourfold degenerate Landau level at Fermi energy is excepted. They 
make density of states exhibit many pairs of asymmetric peak structure and a delta- 
function-like peak, respectively. Among the various factors influencing the density of 
states, only the field strength is important. The strength, but not the period and 
direction, strongly affects the energies of band-edge states, in which the energies of 
local maxima (minima) increase (decrease) in the increase of the field strength. The 
predicted magneto-electronic properties could be examined by measurements on the 
magneto-optical absorption spectra. 
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Figure 1. The primitive unit cell of a monolayer graphene in a perpendicular uniform 
magnetic field and a spatially modulated magnetic field along the armchair direction. 
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Figure 2. The energy bands of (B=40 T, B'=4 T, R=l) at low energy (a) and high 
energy (b) regimes. Those of stronger strength B'=8 T, shorter period R=l/2, longer 
period R=2. along different direction (zigzag), respectively, arc shown in (c), (d), (c) 
and (f). The low energy bands of (B=40 T, B'=0) are also shown in (a) by the dotted 
curves. Notice that b' = \/36 (36) for armchair (zigzag) direction. 
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Figure 3. The low-frequency density of states along armchair direction at (a) (B=40 
T, B'=4 T, R=l), (b) (B=40 T, B'=8 T, R=l), (c) (B=40 T, B'=4 T, R=2), and (d) 
along the zigzag direction at (B=40 T, B'=4 T, R=l). 
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Figure 4. Energies (wc's) of the first six band-edge states, (a) and (b) are their 
dependence on the strength and period respectively. 



